Atom interferometry measurement of the electric polarizability of lithium 
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Using an atom interferometer, we have measured the static electric polarizability of ^Li a = 
(24.33 ± 0.16) X 10"^° m^ = 164.19 ± 1.08 atomic units with a 0.66% uncertainty. Our experiment, 
which is similar to an experiment done on sodium in 1995 by D. Pritchard and co-workers, consists 
in applying an electric field on one of the two interfering beams and measuring the resulting phase- 
shift. With respect to D. Pritchard's experiment, we have made several improvements which are 
described in detail in this paper: the capacitor design is such that the electric field can be calculated 
analytically; the phase sensitivity of our interferometer is substantially better, near 16 mrad/%/Hz; 
finally our interferometer is species selective so that impurities present in our atomic beam (other 
alkali atoms or lithium dimers) do not perturb our measurement. The extreme sensitivity of atom 
interferometry is well illustrated by our experiment: our measurement amounts to measuring a slight 
increase At; of the atom velocity v when it enters the electric field region and our present sensitivity 
is sufficient to detect a variation Ati/w « 6 x 10~^^. 



I. INTRODUCTION 

The measurement of the electric polarizability a of 
an atom is a difficult experiment: this quantity cannot 
be measured by spectroscopy, which can access only to 
polarizability differences, and one should rely either on 
macroscopic quantity measurements such as the electric 
permittivity (or the index of refraction) or on electric 
deflection of an atomic beam. For a review on polariz- 
ability measurements, we refer the reader to the book by 
Kresin and Bonin ij. For alkali atoms, all the accurate 
experiments were based on the deflection of an atomic 
beam by an inhomogeneous electric field and, in the case 
of lithium, the most accurate previous measurement was 
done in 1974 by Bederson and co-workers 0, with the 
following result a — (24.3±0.5) x 10~^° m^. However, 
in 2003, Amini and Gould, using an atomic fountain 
have measured the polarizability of cesium atom with a 
0.14% relative uncertainty, which is presently the small- 
est uncertainty on the electric polarisability of an alkali 
atom. 

Atom interferometry, which can measure any weak 
modification of the atom propagation, is perfectly 
adapted to measure the electric polarizability of an atom: 
this was demonstrated in 1995 by D. Pritchard and co- 
workers with an experiment on sodium atom and they 
obtained a very high accuracy, with a statistical and sys- 
tematic uncertainties both equal to 0.25%. This experi- 
ment was and remains difficult because an electric field 
must be applied on only one of the two interfering beams: 
one must use a capacitor with a thin electrode, a septum, 
which can be inserted between the two atomic beams. 

Using our lithium atom interferometer 0, |^ , we have 
made an experiment very similar to the one of D. 
Pritchard 0] and we have measured the electric polar- 
izability of lithium with a 0.66% uncertainty, limited by 



the uncertainty on the mean atom velocity and not by 
the atom interferometric measurement itself j^. In the 
present paper, we are going to describe in detail our ex- 
periment with emphasis on the improvements with re- 
spect to the experiments of D. Pritchard's group 0, 
we have designed a capacitor with an analytically calcu- 
lable electric field; we have obtained a considerably larger 
phase sensitivity, thanks to a large atomic flux and an ex- 
cellent fringe visibility; finally our interferometer, which 
uses laser diffraction, is species selective: the contribution 
of any impurity (heavier alkali atoms, lithium dimers) to 
the signal can be neglected. 



We may recall that several experiments using atom in- 
terferometers have exhibited a sensitivity to an applied 
electric field [Tol Ull but these experiments were not 
aimed at an accurate measurement of the electric po- 
larizability. Two other atom interferometry experiments 
[13L using an inelastic diffraction process, so that 
the two interfering beams are not in the same internal 
state, have measured the difference of polarizability be- 
tween these two internal states. Finally, two experiments 

S IB 113 
have measured the Aharonov-Casher phase 
: this phase, which results from the application of 
an electric field on an atom with an oriented magnetic 
moment, is proportional to the electric field. 



This paper is organized as follows. We briefly recall the 
principle of the experiment in part II. We then describe 
our electric capacitor in part III and the experiment in 
part IV. The analysis of the experimental data is done in 
part V and we discuss the polarizability result in part VI. 
A conclusion and two appendices complete the paper. 
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FIG. 1: Schematic drawing of our experiment based on 
a Mach-Zehnder atom interferometer; a collimated atomic 
beam, coming from the left, is diffracted by three laser stand- 
ing waves made by reflecting three laser beams on the mirrors 
Ml, M2 and M-j,. The output beam labelled 1 is selected by 
a slit and detected by a hot-wire detector D. The capacitor 
with a septum between the two interfering beams is placed 
just before the second laser standing wave. The x, y and z 
axis are defined. 



II. PRINCIPLE OF THE MEASUREMENT 

If we apply an electric field E on an atom, the energy 
of its ground state decreases by the polarizability term: 
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FIG. 2: Two schematic drawings of the capacitor: top view 
of the capacitor cut in the atomic beam plane (upper panel) 
and end view as seen by an observer located on the atomic 
beam (lower panel). The axis being the same as in figure I, 
the septum and the electrodes are parallel to the y, z plane, 
with the septum at a: = and the electrodes at a; = ±/i. The 
high voltage electrodes at the V = Vo potential extends from 
z — —a to 2 = +a, while the guard electrodes extend outside 
with \z\ > a. The septum and the guard electrodes are at the 
y = potential. 



III. THE ELECTRIC CAPACITOR 
A. Capacitor design 



U = -2TreoaE^ (1) 

When an atom enters a region with a non vanishing 
electric field, its kinetic energy increases by —U and 
its wave vector k becomes k + Afc, with Afc given by 
Afc — 2TTeoaE'^m/ {hk). The resulting phase shift <j) of 
the atomic wave is given by: 

<l> ^ J E\z)dz (2) 

where we have introduced the atom velocity v = hk/m 
and taken into account the spatial dependence of the elec- 
tric field along the atomic path following the z-axis. This 
phase shift is inversely proportional to the atom velocity 
and this dependence will be included in our analysis of 
the results. 

The principle of the experiment, illustrated on figure 
^ is to measure this phase shift by applying an electric 
field on one of the two interfering beams in an atom in- 
terferometer 0. This is possible only if the two beams 
are spatially separated so that a septum can be inserted 
between the two beams. This requirement could be sup- 
pressed by using an electric field with a gradient as in ref- 
erence but it seems difficult to use this arrangement 
for a high accuracy measurement, because an accurate 
knowledge of the values of the field and of its gradient at 
the location of the atomic beams would be needed. 



To make an accurate measurement, we must know pre- 
cisely the electric field along the atomic path and guard 
electrodes are needed so that the length of the capacitor is 
well defined, as discussed by D. Pritchard and co-workers 
PI . It would probably be better to have guard electrodes 
on both electrodes, but it seems very difficult to draw 
guard electrodes on the septum and to put them in place 
very accurately. Therefore, as in reference P|, we have 
guard electrodes only on the massive electrodes. How- 
ever, we have chosen to put our guard electrodes in the 
plane of the high voltage electrode. With this choice, the 
calculation of the electric field can be done analytically. 
Figure |21 presents two schematic drawings of the capac- 
itor and defines our notations, while an artist's view is 
presented in figure Like in reference i, our capacitor 
is as symmetric as possible with respect to the septum 
plane, but, for a given experiment, only one half of the 
capacitor is used, the other part creating no electric field 
with V = everywhere. 



B. Calculation of the capacitor electric field 

1. From three dimensions to two dimensions 

If a dielectric slab with a permittivity > 1 is intro- 
duced in a plane capacitor, the field lines are distorted 
and concentrated towards the slab. Because our capac- 
itor contains dielectric spacers, one could fear a similar 
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FIG. 3: Artist's view of the capacitor: we have also shown 
schematicaUy some wires connecting the electrodes, the im- 
portant point being that these wires are not close to the 
atomic beams. 



effect but this effect does not exist when tlie dielectric 
slab fills completely the gap between the electrodes. Fol- 
lowing figurc[51 the spacers, with a dielectric constant Sr 
completely fill the space \y\ > ?/o, with vacuum in the 
rest of capacitor, \y\ < yo. Let V (respectively Vi) be 
the potential with |y| < yo {\y\ > Vo)- In the |y| = yo 
planes, the continuity equations give: 



V ^Vi and — = 

dy dy 



(3) 



As shown below, a 2D solution of the Laplace equation 
of the form V{x, z) exists for the capacitor. Then clearly, 
if we take Vi = V, this solution fulfills both continuity 
equations on the dielectric borders in the \y\ = ya planes 
and we know that the solution is unique. 



2. Calculation of the potential and of the electric field 

We consider only one half (with x > 0) of the capacitor 
represented on figure |2J We know the potential on the 
borders of the capacitor, with V{x = 0, z) = and V{x — 
h, z) = Vb if \z\ < a and V{x = ft,, z) = if \z\ > a. To 
get the potential everywhere, we start by calculating the 
Fourier transform V{k) of V{x ^ h, z) : 



V{x,z)^ 



1^ 



sinh (kx) 
sinh (fc/i) 



exp (ikz) dk (5) 



from which we can deduce the electric field everywhere. 
On the septum surface x = 0, the electric field is parallel 
to the X-axis and we can calculate exactly the integral 
of E'^ on this surface (see appendix A). We need the 
capacitor effective length L^ff which is defined by: 



Leff - — 



^0 



+ 00 



E^dz 



(6) 



where Eq = Vo/h is the electric field of an infinite plane 
capacitor with the same electrode spacing h. Using equa- 
tion (|2H1), we get the exact value of Leff- 



Leff - 2a 



coth 



/ 7ra \ h 
\ h J ira 



2a (7) 

TT 



where exponentially small corrections of the order of 
exp [—2TTa/h] have been neglected in the approximate re- 
sult. 

However, the atoms do not sample the electric field 
on the septum surface but at a small distance, of the 
order of 50 fim in our experiment, and what we need 
is the integral of E^ along their mean trajectory. This 
mean trajectory is not exactly parallel to the septum, 
but it is easier to calculate the integral along a constant 
X line. We may use either the potential given by equation 
(O or Maxwell's equations to relate the field components 
near the septum surface to its value on the surface. The 
calculation, given in appendix A, proves that the first 
correction to the effective length is proportional to x^: 



L 



eff 



2a - 



2h 2'Kx'^ 



3h 



(8) 



The x^ correction is a fraction of the main term 2a equal 
to 7rx^/(3aft,) and with our dimensions (x ~ 50 /im, /i w 2 
mm and a w 25 mm), this correction is close to 5 x 
10^^ Leff- This correction is negligible at the present 
level of accuracy and we will use the value of L^ff given 
by the approximate form of equation ((TJ . More precisely, 
we will write: 



Vik) 



V2^ J- 



V{x = h, z) exp [—ikz) dz (4) 



Then, for a separable harmonic function of x and z, with 
a z-dependence of the exp [ikz) form, its x-dependence 
is necessarily given by a linear combination of the two 
functions exp (±A:x). Using this result and the conditions 
on the borders at x = and at x = /i, we get the value 
of V{x, z) everywhere: 



E^dz = 



2a 



2 



(9) 



C. Construction of the capacitor 

Let us describe how we build this capacitor. The ex- 
ternal electrodes are made of glass plates (80 mm long in 
the z direction, 35 mm high in the y direction, 10 mm 
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thick in the x direction) covered by an evaporated alu- 
minium layer. To separate the guard electrodes from the 
high voltage electrode, a gap is made in the aluminium 
layer by laser evaporation "ig*!. We found that 100 /im 
wide gaps give a sufficient insulation under vacuum to 
operate the capacitor up to Vb = 500 V. These gaps are 
separated by a distance 2a = 50 mm, so that the two 
guard electrodes are 15 mm long. The glass spacers are 
~ 2 mm thick plates of float glass (10 mm x 80 mm) 
used without further polishing. The distance from the 
spacer inner edge to the atomic beam axis is equal to ~ 7 
mm. 

We have found that a float glass plate is flat within 
±2 /im over the needed surface. This accuracy appeared 
to be sufficient for a first construction, as the main ge- 
ometrical defects are due to the way we assemble the 
various parts by gluing them together and by an imper- 
fect stretching of the septum. We use a double-faced tape 
ARCLAD 7418 (from Adhesive Research) to assemble the 
spacers on the external electrodes. 

The septum is made of 6 /im thick mylar from Goodfel- 
low covered with aluminium on both faces. In a first step, 
the mylar sheet is glued on circular metal support. It is 
then covered by a thin layer of dish soap diluted in water 
and the mylar is heated near 80°C with a hot air gun. 
Then, we clean the mylar surface with water and let it 
dry. After this operation, the mylar is well stretched and 
its surface is very flat. We have measured the resonance 
frequencies of the drum thus formed, from which we de- 
duced a surface tension of the order of 50 N/m (this value 
is only indicative as this experiment was made with an- 
other mylar film which was 20 /im thick). Once stretched, 
the mylar film is glued on one electrode-spacer assembly 
with an epoxy glue EPOTEK 301 (from Epoxy Technol- 
ogy), chosen for its very low viscosity, and then it is cut 
with a scalpel. In a final step, a second electrode-spacer 
assembly is glued on the other face of the mylar. Finally, 
as shown in figure |3| wires are connected to the vari- 
ous electrodes using an electrically conductive adhesive 
EPOTEK EE129-4 (also from Epoxy Technology). 

D. Residual defects of the capacitor 

We are going to discuss the various points by which 
the real capacitor differs from our model. 



1. 2D character of the potential 

We have shown that the potential V is reduced to a 2D 
function when an homogeneous dielectric slab fills com- 
pletely the gap between the electrodes, with the border 
between the vacuum and the slab being a plane perpen- 
dicular to the electrodes. The real dielectric slab is the 
superposition of a tape, a glass spacer and a glue film, 
each material having a different permittivity tr- The 
differences in permittivity perturb the potential which 



should take a 3D character extending on a distance com- 
parable to the tape or glue film thicknesses. This pertur- 
bation seems negligible, because the tape and the glue 
films are very thin and also because these three dielectric 
materials have not very large values. 

2. Do we know the potential everywhere on the border? 

Our calculation assumes that the potential is known 
everywhere on the border. But, on the high-voltage elec- 
trode, we may fear that the potential is not well defined in 
the 100 /im wide dielectric gaps separating the high volt- 
age and guard electrodes as these gaps might get charged 
in an uncontrolled way. This is not likely if the volume 
resistivity p of the pyrex glass used is not too large: more 
precisely, the time constant for the charge equilibration 
on the gap surface is given by e^p (within numerical fac- 
tors of the order of one) and this time constant remains 
below 1 second if p < 10^^ Vl.ia. We have found sev- 
eral values of the resistivity of pyrex glass at ordinary 
temperature, in the range (4 — 8) x 10^ il.m and with 
such a conductivity, this time constant is below 10^^ s. 
Our calculation neglects the surface conductivity, due to 
the adsorbed impurities, which should further reduce this 
time constant. 

Therefore, we think that it is an excellent approxima- 
tion to assume that the potential V makes a smooth 
transition from V = Q to V — Vq in the gaps. Then, 
using equation H2t)|l . it is clear that the detailed shape 
of the transition has no consequence as these details are 
smoothed out by the convolution of V{z) by the function 
g[z) which has a full width at half maximum equal to 
1.12ft,. We can use equation iQ to calculate the effective 
length, provided that we add to the length of the high 
voltage electrode the mean width of the two gaps. In the 
present work, we have taken one gap width, 100 /xm, as 
a conservative error bar on the effective length. A supe- 
riority of our capacitor design is that these gaps are very 
narrow, thus minimizing the corresponding uncertainty 
on the capacitor effective length and we hope to be able 
to further reduce this uncertainty. 



3. Parallelism of the electrodes 

The thickness h of the capacitor, which is the sum of 
the thickness of the spacers, the tape and the glue film, 
is not perfectly constant. Using a Mitutoyo Litematic 
machine, we have measured, with a ±1 /im uncertainty, 
the capacitor thickness as a function of z, in the center 
line of the two spacers, &t y — ±12 mm. The average 
of these two measurements gives the thickness h{z) in 
the y — Q plane around which the atom sample the elec- 
tric field. The thickness h is not perfectly constant but 
it is well represented by a linear function of z, given by 
h{z) = ho + hi(z/a), the maximum deviation hi being 
considerably smaller than the mean h value noted ho- 
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As these deviations are very small (see below), it seems 
reasonable to use use equation Q provided that terms 
involving powers of h are replaced by their correct aver- 
ages. The first term in corresponds to the integral of 
E'^ over the capacitor length, from z — —a to z — a and 
we must take the average value of h^^ over this region. 
Neglecting higher order terms, this average is given by: 



1 



1 + M 

hi 



(10) 



In equation (O, the second term in h~^, corresponds to 
end effects and this quantity must be replaced by the 
following two-point average: 



1 



1 



h{z = —a) h{z = a) 



1 

ho 



hi 



(11) 

Both corrections involve the same factor [1 + (/if//ig)]. 



4- Summary of the capacitor dimensions 

Although the capacitor is as symmetric as possible, 
this symmetry is only approximate and we give the pa- 
rameters for the half we have used for the set of measure- 
ments described below. The length of the high voltage 
electrode, including one gap width is 2a — 50.00 ± 0.10 
mm, the error bar being taken equal to one gap width, 
as discussed above. The distance h between the elec- 
trodes gap width is described by /iq = 2.056 ± 0.003 
mm and hi = 3.2 x 10"'^ mm. The correction term 
hl/hg = 2.4 X 10^^ is completely negligible. 



IV. THE EXPERIMENT 

In this part, we are going to recall the main features 
of our lithium atom interferometer, to give the values of 
various parameters used for the present study, to present 
the data acquisition procedure and the way we extract 
the phases from the data. 



A. Our interferometer 

Our atom interferometer is a Mach-Zehnder interfer- 
ometer using Bragg diffraction on laser standing waves. 
Its design is inspired by the sodium interferometer of D. 
Pritchard and co-workers [2oL I2H and by the metastable 
neon interferometer of Siu Au Lee and co-workers [23| . 
A complete description has been published (SmQ]. 

The lithium atomic beam is a supersonic beam seeded 
in argon and, for the present experiment, we have worked 
with a low argon pressure in the oven pq = 167 mbar, 
because the detected lithium signal increases when the 
argon pressure decreases. The oven body temperature 



is equal to 973 K, fixing the vapor pressure of lithium 
at PLi = 0.8 mbar and the nozzle temperature is equal 
to To = 1073 K. With these source conditions, following 
our detailed analysis 0, the argon and lithium velocity 
distributions are described by a parallel speed ratio for 
argon equal to S^^ ji,, — 8.3 and a parallel speed ratio for 
lithium equal to S\\ i^i = 6.2 (the parallel speed ratio is 
defined by equation (|15|l below). 

We use Bragg diffraction on laser standing waves at 
A = 671 nm: the laser is detuned by about 3 GHz on 
the blue side of the - ^^3/2 transition of the ''Li 

isotope, the signal is almost purely due to this isotope, 
which has a natural abundance equal to 92.41%, and not 
to the other isotope ^Li. Moreover, any other species 
present in the beam, for instance heavier alkali atoms or 
lithium dimers, is not diffracted and does not contribute 
to the signal. 

The case of lithium dimers deserves a special discus- 
sion because they are surely present in the beam and 
the lithium dimer has an absorption band system due its 
A — X transition with many lines around 671 nm: for 
most rovibrational levels of the X state, the absorption 
transition which is closest to the laser frequency has a 
very large detuning of the order of hundreds to thousands 
of GHz and the intensity of this resonance transition 
is also weaker than the resonance transition of lithium 
atom, because of the Franck-Condon factor. Therefore, 
the lithium dimers have a negligible probability of diffrac- 
tion and do not contribute to the interferometer signals. 

The interference signals in a Mach-Zehnder interfer- 
ometer are given by: 



/ = /o[l + VcosV'] 



(12) 



where the phase ijj of the interference fringes can be writ- 
ten: 



■0 = p{2kL){xi -I- ^3 - 2x2) + 4> 



(13) 



The first term of ip is particular to three-grating inter- 
ferometers: the diffracted beam of order p by grating i 
has a phase dependent on the grating position Xi. In 
our case of laser diffraction, the grating position Xi is 
given by the mirror position Mi and 2kL is the grating 
wavevector, where fc^ is the laser wavevector. This phase 
term is very interesting because it is non dispersive and it 
is commonly used to observe interference fringes. In our 
case, we scan the position X3 of mirror M3 by a piezoelec- 
tric translation stage. The second term represent any 
phase difference between the two beams and in particu- 
lar, it will represent the phase shift due to the application 
of an electric field on one of the two paths. 



B. Introduction of the capacitor 

In the present work, we have only used the diffraction 
order p = 1 so that the center of the two beams are 
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FIG. 4: Interference signal observed by displacing the mirror 
Ms: the signal is expressed in counts per second, but the 
counting time is equal to 0.1 second only and the background 
signal has been recorded at the end of the experiment. The 
full curve is the best fit and the fit residuals are plotted, at 
the bottom of the figure. 



separated by about 90 /im at the location of the capaci- 
tor, which is located just before the second laser standing 
wave. The capacitor is attached to the top of the vacuum 
chamber and not to the rail supporting the laser standing 
wave mirrors M^: in this way, we do not increase the vi- 
brations of the mirror positions Xi The capacitor is held 
by a translation stage along the x-direction, which can be 
adjusted manually thanks to a vacuum feedthrough and 
a double stage kinematic mount built in our laboratory. 
The first stage, operated with screws, can be used only 
when the experiment is at atmospheric pressure while the 
second stage, actuated by low-voltage piezo-translators, 
can be adjusted under vacuum. When the septum is 
inserted between the two atomic paths, the atom prop- 
agation is almost not affected by its presence and, as 
shown in figure 01 we have observed a fringe visibility 
equal V = 84 ± 1 % and a negligible reduction of the 
atomic flux. 

To optimize the phase sensitivity (see a discussion in 
reference 0), we have opened the collimation slit 5*1 and 
the detection slit Sd (see reference 6]) with widths ei = 
18 /im and cd — 50 /^m, thus increasing the mean flux up 
to 10*^ counts/s and slightly reducing the fringe visibility 
down to Vo = 62% (see figure Isj. 



C. The data acquisition procedure 

We have made a series of recordings, labelled by an 
index i from 1 to 44, with Vq — when i is odd, and 
with Vq ^ 0, when i is even with Vq « 10 x i Volts. For 
each recording, we apply the same linear ramp on the 
piezo-drive of mirror M3 in order to observe interference 
fringes and 471 data points are recorded with a counting 
time per channel equal to 0.36 s. Figure [3 presents a pair 




50 100 150 200 250 300 350 400 450 500 550 



Channel number 



FIG. 5: Experimental signals corresponding to Vb = (black 
dots) and to Vb ~ 260 Volts (grey dots). Their sinusoidal fits 
using equation 11411 are represented by the full curves. The 
phase shift due to the polarizability effect is close to Stt and 
the fringe visibility is reduced because of the dispersion of the 
phase shift associated to the velocity spread of the lithium 
atoms. 



of consecutive recordings. 

The high voltage power supply has stability close to 
lO"** and the applied voltage is measured by a HP model 
34401A voltmeter with a relative accuracy better than 
10-^ 



D. Extracting phases from the data 

For each recording, the data points Ii{n) have been 
fitted by a function: 



with 'ijji{n) 



loi [1 + Vi cos ^i{n)] 
ai + bin + av? 



(14) 



where n labels the channel number, represents the ini- 
tial phase of the pattern, hi an ideal linear ramp and 
Ci the non-linearity of the piezo-drive. For the Vq = 
recordings, Oj, hi and Ci have been adjusted as well as the 
mean intensity I^i and the visibility Vi. For the V ^ 
recording, we have fitted only a^, I^i and Vi, while fixing 
hi and Ci to their value hi-i and Ci_i from the previous 
V = Q recording. We think that our best phase mea- 
surements are given by the mean phase tpi obtained by 
averaging ipi(n) over the 471 channels. The la error bar 
of these mean phases are of the order of 2 — 3 mrad, in- 
creasing with the applied voltage up to 23 mrad because 
the visibility is considerably lower when the applied volt- 
age Vq is large (see figure El . This rapid decrease of the 
visibility is due to the velocity dependence of the phase 
and to the velocity distribution of the lithium atoms. 
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FIG. 6: Relative fringe visibility V/Vo with Vo = 62% as a 
function of the applied voltage Vo- The points are our mea- 
surements and the full curve is our best fit using equations 



FIG. 7: The mean phase value (tpi), in radian, of the Vb = 
recordings is plotted as a function of the recording starting 
time t in minutes (with an arbitrary origin). The straight 
line is the best linear fit, corresponding to a phase drift of 
7.5 ±0.2 mrad/minute. The fit residuals are also plotted and 
from t — 15 till t = 100 minutes, the residuals exhibit an 
oscillating pattern with a period close to 17 minutes. 



The mean phase values ipi values of the Vb = record- 
ings are plotted in figured they present a drift equal to 
7.5 ± 0.2 mrad/minute and some scatter around this reg- 
ular drift. The most natural explanation for this drift is 
a change of the phase ■0 resulting from a variation of the 
mirror positions Xi : V' changes by 1 radian for a variation 
of {xi + X3 — 1x2) equal to 53 nm. We have verified that 
the observed drift has the right order of magnitude to 
be due to the differential thermal expansion of the struc- 
ture supporting the three mirrors: its temperature was 
steadily drifting at 1.17 x lO^^^ K/minute during the ex- 
periment and the support of mirror M3 differs from the 
other supports, as it includes a piezo translation stage, 
which is replaced by aluminium alloy for the mirrors M\ 
and M2. Presently, we have no explanation of the phase 
scatter, which presents a quasi-periodic structure as a 
function of time: its rms value is equal to 33 milliradian 
and, unfortunately, this scatter gives the dominant con- 
tribution to our phase uncertainty. 

The phase shift ((/)(Vb)) due to the polarizability effect 
(the average () recalls that our experiment makes an av- 
erage over the velocity distribution, as discussed below) 
is taken equal to (0(Vb)) = -0^ — (V^i-i + V'i+i) /2 where 
the recording i corresponds to the applied voltage Vq: the 
average of the mean phase of the two Vq = recordings 
done just before and after is our best estimator of the 
mean phase of the interference signal in zero field. In fig- 
urc|Sl we have plotted the phase shift ^(Vb) as a function 
of the applied voltage Vb . We have chosen for the error 
bar on 0(Vb) the quadratic sum of the Icr error bar given 
by the fit of (V'i) and the 33 milliradians rms deviation 
of the Vb = phase measurements. 



Phase shift (rad) 
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FIG. 8: The measured phase shift (<^(Vo)) in radian is plotted 
as a function of the applied voltage Vq: the best fit using equa- 
tions 14151 is represented by the full curve and the residuals 
are plotted in the lower graph with an expanded scale. 



V. ANALYSIS OF THE SIGNALS: THE EFFECT 
OF THE LITHIUM VELOCITY DISTRIBUTION 

To interpret the experimental data, we must take into 
account the velocity distribution of the lithium atoms. 



A. Velocity averaging of the interference signal 



We assume that this velocity distribution is given by: 



8 



: exp 



{{v-u)Sii/uy 



(15) 



where u is the most probable velocity and S'n is the par- 
allel speed ratio. With respect to the usual form of the 
velocity distribution for supersonic beams, we have omit- 
ted a factor which is traditionally introduced [2^ but 
when the parallel speed ratio S^^ is large enough, this 
has little effect and the main consequence of its omis- 
sion is a slight modification of the values of u and 5|| . 
Here, what we consider is the atoms contributing to the 
interferometer signal and their velocity distribution may 
differ slightly from the velocity of the incident beam, as 
Bragg diffraction is velocity selective. The experimental 
signal can be written: 



Jo jdvPiv) 



I + Vq cos 



/o[l + (V> cos(V.+ (0))] 



(16) 



where is the value of the phase 4> for the velocity 
V = u. If we introduce S = {v — u)/u and expand (u/v) 
in powers of S up to second order, the integral can be 
taken exactly, as discussed in appendix B. However, the 
accuracy of this approximation is not good enough when 
(pm is large and we have used direct numerical integration 
to fit our data. 



B. Numerical fit of the data 

Using equations (|15|) and Ijltill , we have fitted the mea- 
sured phase 4> and visibility V as a function of the applied 
voltage V. The phase measurements received a weight 
inversely proportional to the square of their estimated 
uncertainty and we have adjusted two parameters, the 
value of (pm/V^ and the parallel speed ratio S'||. The re- 
sults of the fits are presented in figures |H1 and |H1 The 
agreement is excellent, in particular for the phase data, 
and we deduce a very accurate value of 4>m/VQ: 



jVf^ = (1.3870 ±0.0010) X 10^"* rad/V^ 



(17) 



where the error bar is equal to la. The relative uncer- 
tainty on (f)m/V^ is very small, 0.072%, which proves the 
quality of our phase measurements. We also get an ac- 
curate determination of the parallel speed ratio: 



S\\ = 8.00 ±0.06 



(18) 



This value of the parallel speed ratio is larger than the 
predicted value for our lithium beam, S'n = 6.2 (see 
above) and this difference can be explained by the veloc- 
ity selective character of Bragg diffraction. 



VI. THE ELECTRIC POLARIZABILITY OF 
LITHIUM 

The lithium electric polarizability is related to the 
value of (jim/V^ by: 



2'KeQa 



2a 2h 

hi TT 



(19) 



All the geometrical parameters a, /iq and hi describing 
the capacitor are known with a good accuracy and we 
deduce a still very accurate value of the ratio of the elec- 
tric polarizability a divided by the mean atom velocity 
u: 



a/u^ (2.283 ±0.008) x 10" 



m^.S 



(20) 



A. Measurement of tlie mean atom velocity 

We have measured the mean atom velocity u by various 
techniques 

We have made measurements by Doppler effect, mea- 
sured either on the laser induced fluorescence signals or 
on the intensity of the atomic beam which is reduced 
by atomic deflection due to photon recoil. In the first 
case, the laser was making an angle close to 49° with the 
atomic beam and it is difhcult to measure this angle with 
sufhcient accuracy. 

In the second case, we have used a laser beam almost 
contra-propagating with the atoms, so that the uncer- 
tainty on the cosine of the angle is negligible. The signal 
appears as an intensity loss on the atomic beam and, the 
loss is not very large because we have used only one laser, 
so that the atoms are rapidly pumped in the other hy- 
perfine state. The experimental signal is shown in figure 
|5] From a fit of this data, we get a value of the mean 
velocity u = 1066.4 ± 8.0 m/s. 

We have also recorded the diffraction probability as a 
function of the Bragg angle, by tilting the mirror forming 
a standing wave. This experiment is similar to the one 
described inour paper (see figure 3) but it is made 
with a lower power density so that only the first order 
diffraction appears. The diffraction is detected by mea- 
suring the intensity of the zero-order atomic beam, as 
shown in figure ^| Using an independent calibration of 
the mirror rotation as a function of the applied voltage on 
the piezo-actuator, we get a measurement of the Bragg 
angle Ob — h/{mu\L) = 79.62 ± 0.63 /irad corresponding 
to u = 1065.0 ±8.4 m/s. 

These two values are very coherent and we may com- 
bine them to give our best estimate of the most probable 
velocity u: 



(1065.7± 5.8)m/s 



(21) 
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FIG. 9: Intensity losses recorded as a function of the laser 
frequency in MHz. The peak A corresponds to a laser beam 
almost contra-propagating with the atoms while the peak cor- 
responds to laser beam perpendicular to the atomic beam. 
The lower curve is the transmission peaks of a Fabry-Perot 
confocal interferometer used for frequency calibration. 



FIG. 10: Plot of the intensity of the atomic beam as a func- 
tion of the voltage applied on the piezoelectric actuator in- 
ducing the rotation angle of mirror M2. When the Bragg 
condition is fulfilled, the direct beam intensity is reduced by 
diffraction in the diffraction orders ±1. The dots represent 
the measured intensity and the curves their best fit. In a sep- 
arate experiment, we have calibrated the rotation angle as a 
function of the piezo voltage, 5.83 ± 0.03 /irad/V. 



We can compare this measurement with the theoreti- 
cal prediction for supersonic expansions. For a pure ar- 
gon beam, in the hniit of an infinite parallel speed ratio, 
theory predicts u — -^/SfesTo/m m/s where Tq is the 
nozzle temperature and m the argon atomic mass. With 
To = 1073 ± 11 K, we get u = 1056.7 ± 5.5 m/s. Three 
small corrections must be made. We must correct for the 
finite value of the argon parallel speed ratio S\\Ar, esti- 
mated to be S'liyif. = 8.3, using the semi-empirical relation 
of Beijerinck and Verster 23] and the associated correc- 
tion reduces the most probable velocity u by a fraction 
equal to 0.75/S'j^^^ ~ 1.09% (this correction is calcu- 
lated in the limit of a vanishing perpendicular tempera- 
ture [12). We must replace the argon atomic mass by a 
weighted mean of the lithium and argon atomic masses 
and, with 0.86 mbar of lithium in 167 mbar of argon, this 
correction increases the velocity by 0.21%. Finally, we 
must take into account the velocity slip effect: the light 
lithium atoms go slightly faster than the argon atoms. 
This difference has been calculated by numerical simu- 
lation by P. A. Skovorodko and this quantity is ex- 
pected to scale like 5"^^^, so that the correction in our 
case is estimated to be 2.42%. We thus predict a most 
probable velocity u = 1073.0 ± 5.6 m/s, where the un- 
certainty comes solely from the uncertainty on the tem- 
perature Tq. This value is in satisfactory agreement with 
our measurements. 



B. The electric polarizability of lithium 

Using the measured value of the most probable velocity 
u, we get the lithium electric polarizability of ^Li: 

a = (24.33 ±0.16) x 10"^" m^ 

= 164.19 ± 1.08 atomic units (22) 

The final uncertainty bar is equal to 0.66%, resulting 
from the quadratic sum of the 0.54% uncertainty on the 
most probable velocity w, the 0.2% uncertainty on the 
effective length of the capacitor, the 2 x 0.15% uncer- 
tainty due to the capacitor spacing ho and the 0.07% 
uncertainty on the interferometric measurement. Unex- 
pectedly, the atom interferometry result has the smallest 
uncertainty! 

Our measurement is a mean of the polarizability of the 
two hyperfine sublevels F — 1 and F = 2 of ''Li. These 
two levels have not exactly the same polarizability. The 
difference Aa = a{F = 2, Mp = 0) - a{F = 1, Mf = 0) 
has been measured with great accuracy by Mowat [2^ . 
We can express this result as a fraction of the mean po- 
larisability, Aa = -(3.0 ± 0.1) x lO^^a. This difference 
is fully negligible with our present accuracy. 

Our result is in excellent agreement with those of the 
previous measurements of a which are considered as be- 
ing reliable. The 1934 measurement of Scheffers and 
Stark 28], which gave a = (12 ± 0.6) x 10"^" m^, is 
generally considered to be incorrect. The first reliable 
measurement is due to Bederson and co-workers j29l] in 
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1961, who obtained a = (20. ± 3.0) x IQ-^o m^, by using 
the E-H gradient balance method. In 1963, Chamberlain 
and Zorn 30] obtained a = (22. ± 2.) x 10"^" m^ by 
measuring the deflection of an atomic beam. Finally, in 
1974, a second experiment was done by Bederson and co- 
workers , using the same E-H gradient balance method 
improved by the calibration on the polarizability of he- 
lium atom in the ^Si metastable state, and they obtained 
the value a = (24.3 ± 0.5) x 10~3° m^. 

We may also compare our measurement with theoret- 
ical results. Many calculations of a have been published 
and one can find a very complete review with 35 quota- 
tions in table 17 of the paper of King [sj published in 
1997. Here is a very brief discussion of the most impor- 
tant results: 

• the first calculation of a was made in 1959 by Dal- 
garno and Kingston [^3 and gave a = 165 a.u.. 
This calculation was rapidly followed by several 
other works. 

• if we forget Hartree-Fock results near 170 a.u., a 
large majority of the published values are in the 
162 — 166 a.u. range. 

• in 1994, Kassimi and Thakkar 33j have made a de- 
tailed study with two important results. They ob- 
tained a fully converged Hartree-Fock value, a — 
169.946 a.u. and this result, far from the exper- 
imental values, proves the importance of electron 
correlation. They also made a series of nth-order 
Moller-Plesset calculations with n = 2, 3 and 4, 
from which they extract their best estimate with 
an error bar, a — 164.2 ± 0.1 a.u.. 

• In 1996, Yan et al. have made an Hylleraas 
calculation, with the final value a = 164.111±0.002 
a.u., this value and its error bar resulting from a 
convergence study. 

Our result is extremely close to these two very accurate 
calculations of a. Two minor effects have not been taken 
into account in these calculations, namely relativistic cor- 
rection and finite nuclear mass correction, but these two 
effects are quite small. 

The relativistic correction on the polarizability has 
been studied by Lim et al. ■ They have made different 
calculations (Hartree-Fock, second order Moller-Plesset, 
coupled cluster CCSD and CCSD(T)) and in all cases, 
their relativistic result is lower than the non-relativistic 
result with a difference in the 0.05 — 0.07 a.u. range. We 
do not quote here their results: even the Hartree-Fock 
value is lower by 0.45 a.u. than the one of reference [33l |. 
because the chosen basis set is too small. 

As far as we know, no calculation of a has been made 
taking into account the finite nuclear mass. An order of 
magnitude of the associated correction should be pven 
by the hydrogenic approximation: a oc [1 -|- (m/M)] and 
in this approximation, the polarisability of ''Li should be 
larger than its Al — oo value by 0.04 a.u. but one should 



not expect this approximation to predict even the sign of 
the correction. A high accuracy calculation of the finite 
mass effect is surely feasible, following the Hylleraas cal- 
culations of Yan and Drake, who have already evaluated 
the finite mass effect on some energies |36jj and oscillator 
strengths of lithium atom. 

VII. CONCLUSION 

We have made a measurement of the electric polariz- 
ability of lithium atom ''Li by atom interferometry and 
we have obtained a = (24.33 ± 0.16) x 10"^° m^ = 
164.19 ± 1.08 a.u. with a 0.66% uncertainty. Our mea- 
surement is in excellent agreement with the most ac- 
curate experimental value obtained by Bederson and 
coworkers ^ in 1974 and we have reduced the uncer- 
tainty by a factor three. Our result is also in excel- 
lent agreement with the best theoretical estimates of this 
quantity due to Kassimi and Thakkar [s^l and to Yan 
et al. |3J|. The neglected corrections (relativistic effect, 
finite nuclear mass effect) should be at least ten times 
smaller than our present error bar. 

Our measurement is the second measurement of an 
electric polarizability by atom interferometry, the pre- 
vious experiment being done on sodium atom by D. 
Pritchard and co-workers in 1995 (see also [2l| and 
0). This long delay is explained by the difficulty of 
running an atom interferometer with spatially separated 
beams. Using a similar experiment, J. P. Toennies and 
co-workers have compared the polarizabilities of helium 
atom and helium dimer but this work is still unpublished 

We want now to insist on the improvements we have 
done with respect to the other measurement of an elec- 
tric polarizability by atom interferometry, due to D. 
Pritchard and co-workers [J : 

• The design of our capacitor permits an analyti- 
cal calculation of the integral along the atomic 
path. This property is important for a better un- 
derstanding of the influence of small geometrical 
defects of the real capacitor. In the present exper- 
iment, the uncertainty on the E'^ integral is equal 
to 0.36% and we think that it is possible to reduce 
this uncertainty near 0.1 % with an improved con- 
struction. 

• We have obtained a very good phase sensitivity of 
our atom interferometer: from our recordings, we 
estimate this phase sensitivity near 16 mrad/VHz. 
The accuracy achieved on phase measurement has 
been limited by the lack of reproducibility of the 
phase between consecutive recordings. We will sta- 
bilize the temperature of the rail supporting the 
three mirrors, hoping thus to improve the phase 
stability. Even if we have not been able to fully 
use our phase sensitivity, we have obtained set of 
phase-shifts measurements exhibiting an excellent 
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consistency and accuracy, as shown by the quaHty 
of the fit of figure|Sland by the accuracy, ±0.072%, 
of the measurement of the quantity 4>„i/Vq . 

• In his thesis [s^, T. D. Roberts reanalyzes the mea- 
surement of the electric polarizability of sodium 
atom made by C. R. Ekstrom et al. Q: he esti- 
mates that a weak contribution of sodium dimers 
to the interference signals can be present as ma- 
terial gratings diffract sodium dimers as well as 
sodium atoms and he estimates that, in the worst 
case, this molecular signal might have introduced 
a systematic error as large as 2% on the sodium 
polarizability result. Our interferometer is species 
selective thanks to the use of laser diffraction and 
this type of error does not exist in our experiment. 
Only lithium atoms are diffracted and even, with 
our choice of laser wavelength, only the ^Li isotope 
contributes to the signal. 

The main limitation on the present measurement of 
the electric polarisability of lithium ''Li comes from the 
uncertainty on the most probable atom velocity u. With 
some exceptions, like the Aharonov-Casher phase shift 
|l8"| which is independent of the atom velocity, a phase- 
shift induced by a perturbation is inversely proportional 
to the atom velocity, at least when a static perturbation is 
applied on one interfering beam. This is a fundamental 
property of atom interferometry and clever techniques 
are needed to overcome this difficulty: 

• T. D. Roberts et al. have developed a way of cor- 
recting the velocity dependence of the phase shift 
by adding another phase shift with an opposite ve- 
locity dependence. They were thus able to observe 
fringes with a good visibility up to very large phase 
shift values. 

• our present results prove that a very accurate mea- 
surement can be made in the presence of an impor- 
tant velocity dispersion without any compensation 
of the associated phase dispersion, but by taking 
into account the velocity distribution in the data 
analysis. 

In these two cases, one must know very accurately a ve- 
locity, the most probable velocity in our case and the ve- 
locity for which the correction phase cancels in the case 
of reference 0. The uncertainty on this velocity may 
finally be the limiting factor for high precision measure- 
ments. Obviously other techniques can be used to solve 
this difRculty. 

Finally, we think that our experiment illustrates well 
two very important properties of atom interferometry: 

• The sensitivity of atom interferometry is a natural 
consequence of the well-known sensitivity of inter- 
ferometry in general, which is further enhanced in 
the case of atom by the extremely small value of the 
de Broglie wavelength. Our phase measurement is 



in fact a direct mesaurement of the increase Aw 
of the atom velocity v when entering the electric 
field. Aw is very simply related to the observed 
phase shift: 



Av 

V 



2tt 



(23) 



This variation is extremely small, with Av/v « 
4 X 10"^ for the largest electric field used in this 
experiment, corresponding to w 25 rad. Our ul- 
timate sensitivity corresponds to a phase w 3 mil- 
liradian which means that we can detect a variation 
Av/v « 6x 10~^^, whereas the velocity distribution 
has a FWHM width equal to 21%! 

• when an atom propagates in the capacitor placed 
in our atom interferometer, its wavefunction sam- 
ples two regions of space separated by a distance 
^ 90 /im with a macroscopic object, the septum, 
lying in between and this situation extends over 75 
microsecond duration, without inducing any loss of 
coherence. This consequence of quantum mechan- 
ics remains very fascinating! 
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IX. APPENDIX A: DETAILED CALCULATION 
OF THE CAPACITOR ELECTRIC FIELD 

The first step is to calculate the Fourier transform V{k) 
of V{x = h, z) defined by equation 10}: 



V{k) 



'2tI J -a 



+a 



Vq exp {—ikz) dz 



2Vo sin (ka) 

X r 



1^ 



(24) 



The potential V[x^ z) given by equation (O can then be 
calculated and, from V{x,z), we can deduce the electric 
field everywhere and in particular, on the septum surface 
where it is parallel to the x-axis: 



E^{x = 0, z) 



72^7- 



+ 00 



V{k)g{k)e^p{ikz)dk (25) 



with g{k) = fc/sinh(fc/i). The electric field Er^{x = 0, z) 
is given by the inverse Fourier transform of the product of 
two functions V{k) and g{k). Therefore, the field Ex{x = 
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0, z) is the convolution of their inverse Fourier transforms 
which are V{x — h, z) and g{z): 



X. APPENDIX B: VELOCITY AVERAGE OF 
THE INTERFERENCE SIGNALS 



E^{x^O,z) 



V{x — h, zi)g{z — zi)dzi (26) 



Using reference Hu\ (equation 6 of paragraph 4.111, p. 
511), we get an cxphcit form of g{z): 



9{z) 



2/i2cosh2(ff) 



(27) 



This result proves that the field decreases asymptotically 
like exp {—n\z\/h), when |z| — a ^ h. We can also get a 
closed form expression of the electric field Ex{x — 0,2:) 
but we may get the integral of E^ without this result, 
simply by using the Parseval-Plancherel theorem and ref- 
erence [i^ (equation 4 of paragraph 3.986, p. 506): 



Ex{x — 0, z)'^dz 



1^ 



coth 



(t) 



(28) 

The atoms sample the electric field at a small distance of 
the septum surface and we need the integral of E"^ along 
their mean trajectory. This mean trajectory is not paral- 
lel to the septum, but it is easier to calculate this integral 
along a constant x line. From Maxwell's equations, one 
gets the first correction terms to the field when x does 
not vanish: 



a;2 Q2 j^^ 

Ex{x,z) « Ex{x = Q,z) - — X -^-^ 



Ez{x,z) « X 



BE, 
dz 



(29) 



where the derivatives are calculated for a; = 0. After an 
integration by parts, one gets: 



-t-CXD 



[El{x,z)+El{x,z)] dz 



-OO 



El{x = 0,z) + 2x'^ 



'dE,- 


2' 


dz 





dz (30) 



where we have kept only the first non vanishing correction 
term in x^ . The calculation of the integral is also done 
with the Parseval-Plancherel theorem and, after some al- 
gebra, we get: 



2x^ 



/• + 00 


'dEx' 


J —OQ 


dz 



dz 



„2 20_ 



Tfcoth(^)-l 



sinh^ {^) 



(31) 



where the approximate result is obtained by neglect- 
ing the exponentially small terms of the order of 
exp [— 27r/i/ (a)]. 



We want to calculate: 

/ = Jo y dvP{v) [l + Vo cos (V + ^rn ^) ] (32) 
with the velocity distribution given by: 



Piv) = 



Sii 



: exp 



{{v-u)Si\/uy 



(33) 



Noting d = (v — u)/u and expanding (u/v) in powers of 
6 up to second order, the integral becomes: 



I/Io = 



Sn 



d6 exp 



X [l + Vocos[ij + (f>m{l-S + S^)]] (34) 
which can be taken exactly: 



I/Io 
(V) 

(0) 



= Vn 



[1 + (V) cos (V' + ((/.))] 

5"!! 



St 



1/4 



exp 



■ arctan 



c2 
■^11 



4 {si + </.2„) 



(35) 
(36) 

(37) 



We have tested this approximation by comparing this 
approximate formula with the result of a computer pro- 
gram, for a parallel speed ratio = 8, corresponding 
to our experimental case. As shown on figure ^2 the 
agreement is very satisfactory, at least with our present 
accuracy on visibility values, with differences of the order 
of 1%. 

We have also studied the difference ((<^) — (j>m) as a 
function of 0,„ and the results are presented in figure 
[T^ This difference can reach large values, for instance 
—0.64 rad when (jjm ~ 25 rad. One may remark that, 
as obvious on equation l|T7|l . the difference {{(j)) — (t)m) is 
linear in (jjm and positive, when (jjm is small, and becomes 
negative and roughly cubic in for larger ipm values. 
If the parallel speed ratio 5*11 is large, S\\ 3> 1, as long as 
the linear term is dominant, the velocity averaged phase 
is given by: 



2^2 



(38) 



and not by (/>,„. With 5*11 = 8, the approximate and 
numerical results are almost equal as long as < 18, 
but their difference increases rapidly for larger 0™ values. 
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being close to 0.05 radian when = 25: even if this 
difference is a very small fraction of (pm, this difference 
is not fully negligible and we have decided not to use the 
approximate analytical results H36|l and H37|) to fit the 
data. 



FIG. 11: Calculation of the average visibility (V) in the case 
= 8 as a function of (prn- The approximate result (equa- 
tion I^SJ is plotted with a dashed line while the result of the 
numerical integration is plotted by the full line. 




FIG. 12: Calculation of the average phase {(j>) in the case 
Sii — 8. The quantity {(j>) — (j)m is plotted as a function of (f>m- 
dashed line for the approximate result (equation 1371 and full 
line for the result of the numerical integration. 
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